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Abstract: A possible explanation for the appearance of light fermions and Higgs bosons 
on the four-dimensional domain wall is proposed. The mechanism of light particle trapping 
is accounted for by a strong self-interaction of five-dimensional pre-quarks. We obtain the 
low-energy effective action which exhibits the invariance under the so called r-symmetry. 
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Then we find a set of vacuum solutions which break that symmetry and the five-dimensional 
translational invariance. One type of those vacuum solutions gives rise to the domain wall 
formation with consequent trapping of light massive fermions and Higgs-like bosons as 
well as massless sterile scalars, the so-called branons. The induced relations between low- 
energy couplings for Yukawa and scalar field interactions allow to make certain predictions 
for light particle masses and couplings themselves, which might provide a signature of 
the higher dimensional origin of particle physics at future experiments. The manifest 
translational symmetry breaking, eventually due to some gravitational and/or matter fields 
in five dimensions, is effectively realized with the help of background scalar defects. As 
a result the branons acquire masses, whereas the ratio of Higgs and fermion (presumably 
top-quark) masses can be reduced towards the values compatible with the present-day 
phenomenology. Since the branons do not couple to fermions and the Higgs bosons do not 
decay into branons, the latter ones are essentially sterile and stable, what makes them the 
natural candidates for the dark matter in the Universe. 
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1. Introduction 



The accommodation of our matter world on a four-dimensional surface - a domain wall 
or a thick 3-brane - in a multi-dimensional space-time with dimension higher than four 
has recently attracted much interest as a theoretical concept promoting novel so- 

lutions to the long standing problems of the Planck mass scale ||, f|], GUT scale |7|. ||, 
SUSY breaking scale |)], [k|, electroweak breaking scale [11]- [14] and fermion mass hierar- 
chy Jl5|, |l~6[| . Respectively, an experimental challenge has been posed for the forthcoming 
collider and non-collider physics programs to discover new particles, such as Kaluza-Klein 
gravitons [18, [ljj, radions and graviscalars [2C, 21], branons p^]-p4]|, sterile neutrinos 



@, ||, ||, |7| etc., together with some other missing energy [^] or missing charge effects 



1 29]. The vast literature on those subjects and their applications is now covered in few 
] . Typically, the thick or fat domain wall formation and the trapping 



review articles [30] 



of low-energy particles on its surface (layer) might be obtained [36|-|35] by a number of par- 
ticular background scalar and/or gravitational fields living in the multi-dimensional bulk - 
see however Ref. |3l^| - the configuration of which does generate zero-energy states localized 
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on the brane. Obviously, the explanation of how such background fields can emerge and in- 
duce the spontaneous symmetry breaking is to be found and the domain wall creation, due 
to self-interaction of certain particles in the bulk with low-energy counterparts populating 
our world, may be one conceivable possibility. 

In this paper we consider and explore the non-compact 4 + ^-dimensional fermion 
model with strong local four-fermion interaction that leads to the discrete symmetry break- 
ing, owing to which the domain wall pattern of the vacuum state just arises and allows the 
light massive Dirac particles to live in four dimensions 1 . In such a model the scalar fields 



appear to be as composite ones. We are aware of the possible important role [|3l|, |37|] of a 
non-trivial gravitational background, of propagating gravitons and gauge fields for issues 
of stability of the domain wall induced by a fermion self-interaction. Nonetheless, in order 
to keep track of the main dynamical mechanism, we simplify herein the fermion model just 
neglecting all gravitational and gauge field interaction and, in this sense, our model might 
be thought as a sector of the full Domain- Wall Standard Model. However, we believe that 
the simplified model we treat in the present investigation will be able to give us the plenty 
of quantitative hints on the relationships among physical characteristics of light particles 
trapped on the brane. 

Let us elucidate the domain wall phenomenology in more details and start from the 
model of one five-dimensional fermion bi-spinor ip{X) coupled to a scalar field §{X). The 
extra-dimension coordinate is assumed to be space-like, 

(X a ) = (xp, z) , (Xfj) = (x ,x 1 ,x 2 ,x 3 ) , (g aa ) = (+,-,-,-,-) 

and the subspace of eventually corresponds to the four-dimensional Minkowski space. 
The extra-dimension size is supposed to be infinite (or large enough). The fermion wave 
function is then described by the Dirac equation 

[i la d a - <5>{X)]^{X) =0 , 7a = (7M»-*75), {7«, 7/?} =^9aP , (1-1) 

with 7 Q being a set of four-dimensional Dirac matrices in the chiral representation. 

The trapping of light fermions on a four-dimensional hyper-plane - the domain wall 
- localized in the fifth dimension at z = zq can be promoted by a certain topological, 
z-dependent background configuration of the scalar field (&(X))o = (p(z), due to the ap- 
pearance of zero-modes in the four-dimensional fermion spectrum. In this case, from the 
viewpoint of four-dimensional space-time, Eq. ( |1.1[ ) just characterizes the infinite set of 
fermions with different masses that it is easier to see after the Fock-Schwinger transforma- 
tion, 

[i la d a + <p{z)][i la &* - <p(z) MX) = (-8^ - m 2 z )^{X) ; 

m\ = -d\ + v\z) - 75 ^(z) = m\P L + m 2 _P R , (1.2) 



x One can find some relationship of this mechanism of domain- wall generation to that one of the Top- 
Mode Standard Model E(J used to supply the top-quark with a large mass in four dimensions due to 
quark condensation uniform in the space-time. However, in our case, the vacuum state will receive a scalar 
background defect breaking translational invariance. As well, our model is taken five-dimensional and non- 
compact as compared to six- or eight-dimensional extensions of the Top-Mode Standard Model ^ with 
compact extra dimensions and an essential role played by Kaluza-Klein gravitons and/or gauge fields. 
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where Pl.r = g(l ±75) ar e projectors on the left- and right-handed states. Thus the mass 
operator consists of two chiral superpartners - in the sense of supersymmetric quantum 
mechanics [43]- [45] 

ml = -d 2 + <p 2 (z) T <f/(z) = [-d z ± <p(z) }[d z ± <p(z) ] ; (1.3) 
m\ q + = q + fh 2 ^, fh 2 ^ q~ = q~ m\ , g ± = =fd z + ip(z) ■ (1-4) 



The factorization ( |1.3j ) guarantees the positivity of the mass operator - i.e. the absence 
of tachyons - and the supersymmetry realized by the intertwining relations fll.4| ) entails 
the equivalence of the spectra between different chiralities for non-vanishing masses. As a 
consequence, the related left- and right-handed spinors can be assembled into the bi-spinor 
describing a massive Dirac particle which, however, is not necessarily localized at any point 
of the extra-dimension if the field configuration (p(z) is asymptotically constant. Indeed the 
massive states will typically belong to the continuous spectrum - or to a Kaluza-Klein tower 
for the compact fifth dimension - and spread out the whole extra-dimension. Meanwhile, 
the spectral equivalence may be broken just by one single state, i.e. a proper normalizable 
zero mode of one of the mass operators m\. Let us assume to get it in the spectrum of 
fh 2 + : then from Eq.s (|L~3|) and (|1.4|) it follows that 



q %(x,z) = 0, %(x,z)=^ L (x)ex.p^-J dwip(w)j , (1.5) 

where = Pii^ix) is a free Weyl spinor in the four-dimensional Minkowski space. 

Evidently, if a scalar field configuration has the following asymptotic behavior: namely, 

tp{z) ±C±\z\ v± , Rev± > -1 , C±>0, 

then the wave function ip^ (x, z) is normalizable on the z axis and the corresponding left- 
handed fermion is a massless Weyl particle localized in the vicinity of a four-dimensional 
domain wall. It is easy to check that in this case the superpartner mass operator m 2 _ does 
not possess a massless normalizable solution and if (f(z) is asymptotically constant, with 
C± > and v± = 0, there is always a gap for the massive Dirac states. Further on we 
restrict ourselves to this scenario. 

The important example of such a topological configuration can be derived for the 
system having the free mass spectrum - continuum or Kaluza-Klein tower - for one of 
the chiralities, say, for right-handed fermions. This is realized by a kink-like scalar field 
background 

ip + = Mtanh(Mz) . (1.6) 
The two mass operators have the following potentials 



d 2 z + M 2 [ 1 - 2sech 2 (Mz) ] ; m 2 _ = -d 2 + M 2 , (1.7) 



and the left-handed normalized zero-mode is properly localized around z = 0, in such a 
way that we can set 



V^(x, z) = if> L (x) Mz) , Mz) = W sech(Mz) . (1.8) 
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As a consequence, the threshold for the continuum is at M 2 and the heavy Dirac particles 
may have any masses m > M, the corresponding wave-functions being completely de- 
localized in the extra-dimension. 

On the one hand, if we investigate the scattering and decay of particles with energies 
considerably smaller than M, we never reach the threshold of creation of heavy Dirac 
particles with m > M and all physics interplays on the four-dimensional domain wall with 
thickness ~ 1/M. On the other hand, at extremely high energies, certain heavy fermions 
will appear in and disappear from our world. 

It turns out that the real fermions - quarks and leptons living on the domain wall by 
assumption - are mainly massive. Therefore, for each light fermion on a brane one needs 
at least two five-dimensional proto- fermions ipi(X), ip2(X) in order to generate left- and 
right-handed parts of a four-dimensional Dirac bi-spinor as zero modes. Those fermions 
have clearly to couple with opposite charges to the scalar field §(X), in order to produce 
the required zero modes with different chiralities when {<&(X))q = c/3 + (z): namely, 

[ip-T 3 $(x)]*(x) = , @ = %d a , = 

where 7 a = ^ a <S> I2 are Dirac matrices and r a = I4 (g> a a , a = 1, 2, 3 are the generalizations 
of the Pauli matrices a a acting on the bi-spinor components ipi(X). 

In this way one obtains a massless Dirac particle on the brane and the next task is to 
supply it with a light mass. As the mass operator mixes left- and right-handed components 
of the four-dimensional fermion it is embedded in the Dirac operator (|1.9|) with the mixing 
matrix T\m,f of the fields ipi{X) and ^2 (A). At last, following the general Standard Model 
mechanism of fermion mass generation by means of the Higgs scalars, one can introduce 
the second scalar field H(x) to make this job, replacing the bare mass T\m,f — * t\H(x). 
Both scalar fields might be dynamical indeed and their self-interaction should justify the 
spontaneous symmetry breaking by certain classical configurations allocating light massive 
fermions on the domain wall. From the previous discussion it follows that the minimal set of 
five-dimensional fermions has to include two Dirac fermions coupled to scalar backgrounds 
of opposite signs. In addition to the trapping scalar field, another one is in order to supply 
light domain wall fermions with a mass. Thus we introduce two types of four-fermion 
self-interactions to reveal two composite scalar fields with a proper coupling to fermions. 
As we shall see, these two scalar fields will acquire mass spectra similar to fermions with 
light counterparts located on the domain wall. The dynamical mechanism of creation of 
domain wall particles turns out to be quite economical and few predictions on masses and 
decay constants of fermion and boson particles will be derived. 

2. Fermion model with self-interaction in 5D 

Let us consider the model with the following Lagrange density 

£( 5 W,*) =V + ^ (*r 3 *) 2 + (*n*) 2 , (2.1) 



MX) 
I MX) 



(1.9) 
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where ^f(X) is an eight-component five-dimensional fermion field, see Eq. fll.SD - either 
a bi-spinor in a four-dimensional theory or a spinor in a six-dimensional theory - which 
may also realize a flavor and color multiplet with the total number N = Nf N c of states. 
The ultraviolet cut-off scale A bounds fermion momenta, as the four-fermion interaction 
is supposed here to be an effective one, whereas g\ and g 2 are suitable dimensionless and 
eventually scale dependent effective couplings. 

This Lagrange density can be more transparently treated with the help of a pair of 
auxiliary scalar fields &(X) and H(X), which eventually will allow to trap a light fermion 
on the domain wall and to supply it with a mass 

£( 5 ) If) = 

/VA 3 ATA 3 

m - r 3 $ - nif)* - — $ 2 - — H 2 . (2.2) 

9i 92 

In this model the invariance (when it holds) under discrete r-symmetry transformations 

q, — > Tl ^ ; $ — ► (2.3) 
$ — >t 2 ^; <!>,H — ► (2.4) 
#^T 3 tf; H—>-H, (2.5) 

does not allow the fermions to acquire a mass and prevents a breaking of translational 
invariance in the perturbation theory. This r-symmetry can be associated to the chiral 
symmetry in four dimensions 2 . 

However for sufficiently strong couplings, this system undergoes the phase transition 
to the state in which the condensation of fermion-anti-fermion pairs does spontaneously 
break - partially or completely - the r-symmetry. 

The physical meaning of the scale A is that of a compositeness scale for heavy scalar 
bosons emerging after the breakdown of the r-symmetry. In order to develop the infrared 
phenomenon of r-symmetry breaking, the effective Lagrange density containing the essen- 
tial low-energy degrees of freedom has to be derived. 

To this concern we proceed - only in this Section - to the transition to the Euclidean 
space, where the invariant four-momentum cut-off can be unambiguously implemented. 
Within this framework, the notion of low-energy is referred to momenta \p\ < Aq as com- 
pared to the cut-off A 3> Ao- However, after the elaboration of the domain wall vacuum, 
we will search for the fermion states with masses rrif much lighter than the dynamic scale 
Ao, i.e. for the ultralow-energy physics. Thus, eventually, there are three scales in the 
present model in order to implement the domain wall particle trapping. 

Let us decompose the momentum space fermion fields into their high-energy part 
Whip) = ^(p)$(\p\ - Ao)#(A - \p\), their low-energy part ^i(p) = ^>(p)fl(A - \p\) and 
integrate out the high-energy part of the fermion spectrum, $(£) being the usual Heaviside 
step distribution. More rigorously, the above decomposition of the fermion spectrum should 

2 It can be also related to the chiral symmetry in a six-dimensional space-time where from our five- 
dimensional model can be derived by dimensional reduction. 



- 5 - 



be done covariantly, i.e. in terms of the full Euclidean Dirac operator, 

D = i(@ + r 3 <f> + nH) . (2.6) 

Nevertheless, as we want to concentrate ourselves on the triggering of r-symmetry break- 
ing by fermion condensation, we can safely assume to neglect further on the high-energy 
components of the auxiliary boson fields, what is equivalent to the use of the mean-field or 
large N approximations. Then the low-energy Lagrange density, which solely accounts for 
low-energy fermion and boson fields, can be written as the sum of the expression in Eq. 
( |2.2| ) and the one-loop effective Lagrange density: namely, 

£®(*l, $, H) = 4 5) (^, H) + A£( 5 )(<D, H), (2.7) 

where the tree-level low-energy Euclidean Lagrange density reads 

4 5) =M l (0 + T 3 <S> + T 1 H)y l + — <S> 2 + —H 2 . (2.8) 

9i 92 

The one-loop contribution of high-energy fermions is given by 

A£^(<S>, H) = C - (N/2) ti(X\A\X) , (2.9) 
A = tf(A 2 - D*D) In - (A 2 , - D ] D) In , 

where the normalization constant C is such that AC^ (0, 0) = and the tr operation 
stands for the trace over spinor and internal degrees of freedom. In Eq. fl2.9p the choice of 
normalizations in the logarithms ensures the continuity of spectral density at the positions 
of cut-offs A and Ao- Thereby the spectral flow through the spectral boundary is suppressed. 
In the latter operator A we have incorporated the cut-offs which select out the high-energy 
region defined above [j46[| . 

From the conjugation property D* = T2DT2, it follows that the Euclidean Dirac operator 
is a normal operator, which has to be implemented in order to get a real effective action 
and to define the spectral cut-offs with the help of the positive operator 

D*D = -d^ - d 2 z + M 2 (X) = -d a d a + M 2 (X) , 
M 2 (X) = $ 2 (X) + H 2 (X) - r 3 $$(X) - n 0H(X) . (2.10) 

One can see that, in fact, the scale anomaly only contributes into AC^ , i.e. that part 
which depends on the scales. Thus, equivalently, 

A£^($, H) = C + N [ A ^ tr{X\${Q 2 - D*D)\X) . (2.11) 

A Q 

As we assume that the scalar fields carry momenta much smaller than both scales A S> Aq , 
then the diagonal matrix element in the RHS of Eq. ( 2.11 ) can be calculated either with 
the help of the derivative expansion of the master representation [46] 

.2 n tnMv\ _ ™ [ + °° dt exp{it} 



-00 2?" 



2ni t — ie 



/d n k ( it 1 

^— exp |- ^ [ k 2 + 2iQk a d a - d a d a + M 2 (X) ] I , (2.12) 
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where n is a number of dimensions of the Euclidean space, or by means of the heat kernel 
asymptotic expansion. For n < 5, only three heat kernel coefficients at most are propor- 
tional to non-negative powers of the large parameter Q (see Appendix A) 



(X\$(Q 2 - D*D)\X) = C Q n + dQ n - 2 M\X) 
+Q n ~ A {C 2 [M 2 (X)f + C 3 d a d a [M 2 (X)]} + O (Q n ~ 6 ) , 



(2.13) 



where, for n < 6 and large scales Aq < Q < A, the neglected terms rapidly vanish. For the 
given operator "d(Q 2 — D*D), the coefficients C\ take the following values 



Co 
Ci 
C 2 
C 3 



n 2 n ~ 1 W 2 T(n/2) ' 

1 n 
~2 
n(n 



2"W 2 r(n/2) 
n- 2 

2 n+2 rfn/2Y(n/2) ~ 

n-2 

~3-2 n + 2 n n / 2 T(n/2) 



C <0 ; 

— C > ; 
n(n — 2" 



24 



C <0 , 



(2.14) 



where T(y) stands for the Euler's gamma function. For different possible definitions of 
the effective Lagrange density, involving operators and regularization functions other than 
those ones of Eq. ( 2.11 ), one might obtain in general different sets of Q, albeit their signs 
are definitely firm. As we shall see further on, the negative sign of C\ catalyzes the chiral 
symmetry breaking at sufficiently strong coupling constants. Taking the little trace and 
integrating the RHS of Eq. ( 2.13| ), from Eq. ( [2.11|) one finds, up to a total n-divergence and 
setting C = C (A% - A n )N/n, 



where 



Ax 
A 2 



A ~°° - Ax (A n - 2 -A^~ 2 ) (& + H 2 ) 

+ A 2 (A"" 4 - A™" 4 ) [ {d a <S>f + (d a H) 2 
+ A 2 (A n " 4 - A™" 4 ) [ ($ 2 + H 2 f~ 



N n|5 N 

(n-2)2"" 3 W 2 r(n/2) * 97T 3 ' 

N(n - 2) «T5 N 
(n-4)2 n - 1 W 2 r(n/2) * I^ 3 ' 



(2.15) 

(2.16) 
(2.17) 



As Ao <C A one can neglect the Ao-dependence in Eq. ( 2.15 ) for n > 4, whereas near four 
dimensions the pole in A 2 together with the cut-off dependent factor generate the coefficient 
~ In (A/Aq). For n = 5 we eventually find 



A— >oo 



AA 3 
9vr 3 



;d, 2 + jf / 2 ) + ^(cD 2 + F 2 ) 2 . (2.18) 
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Although the actual values of the coefficients A% and A 2 might be regulator-dependent, as 
already noticed, the coefficients of the kinetic and quartic terms of the effective Lagrange 
density are definitely equal, no matter how the latter is obtained from the basic Dirac 
operator of Eq. ( |2.6| ), This very fact is at the origin of the famous Nambu relation between 
the dynamical mass of a fermion and the mass of a scalar bound state in the regime of 
r-symmetry breaking, as we will see in the next Section. 



3. r-symmetry breaking 

The interplay between different operators in the low-energy Lagrange density (|2.7| ) may 
lead to two different dynamical regimes, depending on whether the r-symmetry is broken 
or not. Indeed, going back to the Minkowski five-dimensional space-time, the low-energy 
Lagrange density can be suitably cast in the form 

£®(*l, <&, H) = %(X)[i? - r 3 $(X) - Tl H{X) ]%{X) 

+ ^d a $(X)d a $(X) + ^d a H(X)d a H(X) 

- ^[<$>\ X )+H 2 (X)) 2 . (3.1) 

We notice that the quartic operator in the potential is symmetric and positive, whilst the 
quadratic operators have different signs for weak and strong couplings, the critical values 
for both couplings being the same within the finite-mode cut-off regularization of Eq. (|2.1§| ) : 
namely, 

g? = 9Tr 3 , i = l,2. (3.2) 

Let us introduce two mass scales Aj in order to parameterize the deviations from the critical 
point 

2A 2 

A l {g l ) = — (g l -gf) . (3.3) 
Taking into account that gi > we have 

Ai( gi ) < (2/9)A 2 (3.4) 
and the effective Lagrange density for the scalar fields takes the simplified form 

4cL(*> H ) = ^ [( d ^f + {d a Hf + 2A!* 2 + 2A 2 tf 2 - (<D 2 + H 2 ) 2 ] . (3.5) 



So far two constants gi, and thereby Aj, play an equivalent role and the related vertices 
are invariant under the replacements *$>i * — * T2^i and subsequent reflection of the scalar 
fields - see Eq. ( |2.5| ). Therefore, without loss of generality, one can always choose 

Ai( 5 i) > A 2 (g 2 ) . (3.6) 
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Whenever both couplings gi are within the range < gi < gf = 97r 3 , then we have 
Ai(gi) < and consequently the potential has a unique symmetric minimum. If instead 
one at least of the couplings gi does exceed its critical value, then the symmetric extremum 
at <3? = H = is no longer a minimum, though either a saddle point for £±2(92) < < Ai (gi) 
or even a maximum for < £.2(92) < ^1(91) ■ If Ai(<?i) > 0, then the true minima appear 
at a non- vanishing vacuum expectation value of the scalar field &(X): namely, 



(/) * 7 = ($(X)) = ±7a7(^) , Hj = (H(X)) 

This follows from the stationary point conditions for constant fields 

[Ai(<7i) - <I> 2 - H 2 ] $ = , 



. 



(3.7) 



[A 2 (<72 



$ 2 1 H 







(3.8) 



and from the positive definiteness of the second variation of the boson effective action for 
constant boson fields. As a matter of fact, if we set 



S 

V[S] 

we readily find 

V[S] - V [S T ] : 



(Si, 5 2 

NA 
4tt 3 



Sj= ±VAi,0 



j d 5 X {-2AiSf(X) + [S t (X)S t (X)} 2 } , 



(3.9) 



ljd 5 xj d 5 Y Si (X) Sj (Y) 



5 2 V[S] 



NA 

4-7T 3 



SSi(X)SSj(Y) 
J d 5 X Si (X) Sj (X) [M/j . + ... , 



+ 



5=5/ 



where 



whereas the Hessian mass matrix M/: 



4Ai 
2Ai - 





1 


i = 1 


2, 






' M 2 


' 


2A 2 ( 










(3.10) 
(3.11) 

(3.12) 



is manifestly positive definite and determines the mass spectrum of the five-dimensional 
scalar excitations. 

A further constant solution of Eq. ( |3.8[ ) does exist for A 2 > 0, i.e. (&)o = and 
(H)o = ±y/A2- However, it corresponds to a saddle point of the potential, as it can be 
seen from Eq. ( |3.10| ) for A\ > A 2 . Likewise, if Ai > A 2 > 0, then the matrix M is 
negative definite at the symmetric point ( < l ) )o = (H)o = which corresponds thereby to a 
maximum. The degenerate situation - i.e. the valley - actually occurs for Ai = A 2 > 0, 
when the rotational r 2 -symmetry is achieved by the Lagrange density but is spontaneously 
broken. The massless scalar state in Eq. ( |3.12| ) arises in full accordance with the Goldstone's 
theorem. 
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The corresponding dynamical effect for the fermion model of Eq. ( |2,l| ) consists in the 
formation of a fermion condensate and the generation of a dynamical fermion mass M - 
see Eq.s (|2.2| ) and ( |2,71 ) - that breaks the T2-symmetry. Its ratio to the heaviest scalar 
mass just obeys the Nambu relation 

Mi = 2M > M 2 , M = (*) = , (3-13) 

the second, lighter composite scalar being a pseudo-Goldstone state. We notice that the 
above relationship holds true independently of the specific values of the coefficients Ai in 
Eq. ( p. 15 ) and, consequently, it takes place in four and five dimensions. However, if we 



properly re-scale the scalar fields according to 

$(X)V£ = ±Ms/A + vu(X) , H{X)^ = vv(X) , v = y / 2ir 3 /N , (3.14) 

in such a way that dim [it] = dim[u] = 3/2, then the low-energy Lagrange density ( |3.1[) can 
be suitably recast in the form 

+ {l/2)d a u(X)d a u(X) - 2M 2 u 2 (X) 

+ {l/2)d a v(X)d a v(X) - 2{M 2 - A 2 )v 2 (X) 

- uA- 1 / 2 [%{X)^ t {X)u(X) -%(X) n %{X)v(X)] 
T 2vMtC x l 2 [u 2 {X) + v 2 {X)]u{X) 

- {v 2 /2)K- 1 [u 2 {X) +v 2 {X)} 2 . (3.15) 

As a consequence, one can see that all the low-energy effective couplings for fermion and 
boson fields do rapidly vanish in the large cut-off limit. In particular, if the cut-off A is 
much larger than the energy range of our physics, then we are dealing with a theory of 
practically free, non-interacting particles. Moreover, this pattern of r-symmetry breaking 
does not provide the desired trapping on a domain wall: the heavy fermions and bosons live 
essentially in the whole five-dimensional space. From now on we shall proceed to consider 
another type of vacuum solutions, which break the five-dimensional translational invariance 
and give rise to the formation of domain walls. 

4. Domain walls: massless phase 

The existence of two minima in the potential of Eq. gives rise to another set of vacuum 



solutions |47]] - [51], which connect smoothly the minima owing to the kink-like shape of 
Eq. ([O]) with M = \/Ai. On variational and geometrical grounds one could expect that 
certain minimal solutions are collinear, just breaking the translational invariance in one 
direction. We specify this direction along the fifth coordinate z. Then one can discover 
two types of competitive solutions p9[-|5l]], 

(J) ($(A)) = ®j{z) = ±Mtanh(Mz) , 

(H(X)) = Hj(z) = ; (4.1) 

(JO ($(A)) = &k{z) = ±Mtanh(/3z) , 

(H(X)) = H K (z) = ±/isech(/3z) . (4.2) 
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Further on we select out only positive signs in vacuum configurations to analyze the scalar 
fluctuations around them, having in mind that our analysis is absolutely identical around 
other configurations. When we insert the second solution {K) into the stationary point 
conditions 

2 [ M 2 - $ 2 - H 2 } $ = d a d a $ , 

2 [ A 2 - H 2 - $ 2 ] H = d a d a H , (4.3) 

we find 

H = y / 2A 2 - M 2 , 13 = VM 2 - fi 2 . (4.4) 

The solution (K) exists only for A 2 < M 2 < 2A 2 and it coincides with the extremum (J) 
in the limit A 2 — ► M 2 /2, /j, — ► 0, j3 — > M. The question arises about which one of the two 
solutions is a true minimum and whether they could coexist if fj, > 0. The answer can be 
obtained from the analysis of the second variation of the bosonic low-energy effective action. 
The corresponding relevant second order differential operator can be suitably written in 
terms of the notations introduced in Eq. (|3.9|): namely, 



[D^ ) = - 5 lj U x - \m[S{X)] ) , U x = , (4.5) 

[M[S(X)] ) = 5ij [-d 2 - 2Ai + 2S k (X)S k (X)} + ASi(X)Sj(X) . (4.6) 

The stationary point solutions (J) and (K) are true minima iff the matrix-valued mass 
operator M[S'(X)] becomes positive semi-definite at the extrema 

Sj = (MtanhMz, 0) , S K = (Mtanh/3z, \x sech(3z) . (4.7) 

Now, at the stationary point (J) the matrix-valued mass operator 

Mj(z) = M[Sj(z)} (4.8) 

turns out to be diagonal with entries 

M J = ~ d z + 4Af2 - 6Af2 sech 2 (Mz) , (4.9) 
Mf = -d 2 + 2M 2 - 2A 2 - 2M 2 sech 2 (Mz) , (4.10) 
M^ 2 = M 2 / = 0. (4.11) 

Both components do represent one-dimensional Schrodinger-like operators, the eigenvalue 
problem of which can be exactly solved analytically. The Schrodinger-like operators ( |4.S| ) 



and ( 4.10 ) can be presented in the factorized form similar to that one of Eq. (|1.3|): namely, 

M^ 1 = \-d z + 2Mtanh(Mz)] [d z + 2Mtanh(Mz)] ; (4.12) 

M^ 2 = M 2 — 2A 2 + [-d z + Mtanh(Mz)] [d z + Mtanh(Mz)] (4.13) 

Therefrom it is straightforward to check that the ground states of the operator Mj(z) are 
described by the real, node-less in z and normalized wave functions 

My <t>j(z) = , M 22 hj(z) = m 2 h hj(z) , m 2 =M 2 -2A 2 ; (4.14) 

4>j(z) = sech 2 (Mz)y / 3M/4 , hj(z) = sech(Mz)^M/2 ; (4.15) 
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in such a way that we can suitably parameterize the shifts of the scalar field with respect 
to the background vacuum solution (J) - see Eq.s ( |3.11| ) and ( f4,l|) - by the following two 
eigenstates of the mass matrix Mj(z): namely, 







A 



n h j(x) = h(x) 



o 

hj{z) 



A 



(4.16) 



where (j)(x), h(x) do eventually represent the ultralow-energy scalar fields on the Minkowski 
space-time, as we shall better see below on. As a consequence, the spectrum of the second 
variation is positive if M 2 > 2A 2 and, in this case, the solution (K) does not exist whilst 
the scalar lightest states are localized on the domain wall. More precisely - see Appendix 
B - the first boson <I> has two states on the brane: a massless state and a heavy massive 
state of mass y/3M. The existence of a massless scalar state around the kink configuration 
(J) is a consequence of the spontaneous breaking of translational invariance - see next 
Section. Other heavy states belong to the continuous part of the spectrum with threshold 
at 2M. The second boson H has only one state on the brane of mass VM 2 - 2A 2 and its 
continuous spectrum starts at \/2M 2 — 2A 2 > M. 

Since the vacuum expectation value of the scalar field (3>(X))o = Mtanh(Mz) has a 
kink shape, its coupling to fermions induces the trapping of the lightest, massless fermion 
state on the domain wall: namely, 



iPil(x] 



ifj (z) , Mz) = sech(Mz) v / M/2 , 



(4.17) 



see Eq.s (|l]^) and ( |1.9| ). The continuum of the heavy fermion states begins at M and 
involves pairs of heavy Dirac fermions. 

In conclusion, at ultralow energies much smaller than M, the physics in the neighbor- 
hood of the vacuum (J) is essentially four-dimensional in the fermion and boson sectors. 
It is described by the massless Dirac fermion 



ip(x) 



(4.18) 



with and ipn(x) being the two-component non-trivial parts of iPil(x) and iI>2r(x) 

respectively, in such a way that we can set 



MX) 



i/j (z) . 



(4.19) 



ip L (x) 
y iPr(x) t 

In result one has two four-dimensional scalar bosons, a massless one and a massive one, 
provided M 2 — 2A 2 <C M 2 otherwise there is decoupling. 

The matrix T3 does not mix the two types of fermions ipiL and ip2R, but the related 
Yukawa vertex in the Lagrange density ( |2.2p mixes left- and right-handed components of 
each of them. As a consequence the massless scalar field 4>(x) does not couple directly to a 
light fermion-anti-fermion pair. Its coupling to fermions involves inevitably heavy fermion 
degrees of freedom. Therefore, the ultralow-energy effective action does not contain a 
Yukawa-type vertex for the field 4>(x), which appears thereby to be sterile. 
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On the other hand, the interaction between light fermions and the second scalar field 
h(x) on the domain wall does achieve the conventional Yukawa form. Indeed, once they 
are projected on the zero-mode space of ijj(x)ij)o(z), the matrices n, ti and T3 act as the 
Dirac matrices in the chiral representation: namely, 



n — > 70 = L J ; t 3 — > 75 = _ . ; r 2 — ► 47075 ; (4.20) 





where the 2x2 unit matrix 1 acts on Weyl components iPl(x) and ipn(x). Meanwhile the 
matrices in the kinetic part of the Dirac operator are projected onto the zero-mode space 
as 

7o% — ► ( ^ ° J = 7o7fc , (4.21) 
y -a k J 

where the Pauli matrices act on the two-component spinors iI>l(%) and iPr(x). As a 
consequence, the matrix t\ — ► 70 just induces the Yukawa mass-like vertex in the effective 
action on the domain wall. 

As a final result, in the vicinity of the vacuum solution (J) of Eq. (|4.l| ), the ultralow- 
energy effective Lagrange density for the light states on the four-dimensional Minkowski 
space-time comes out from Eq.s 0, (|415l) , (p^) , (p^|) and reads 3 



dzLfl ^(I),^^),^^),^^) (4.22) 

-00 



+ l -d^h{x)d»h{x) - l -m 2 h h 2 {x) 



— gfip(x)ip(x)h(x) — Ai(/> (x) — \2(/) (x)h (x) — A3/1 (x) , 

with the scalar mass m| = (M 2 — 2A2) and the ultra- low energy effective couplings given 
by 



2?r 3 f + °° 0t , vr /Mvr 3 



9f = \ TfI dz h j( z )M z ) 



A^A./_oo 4 V AN 



\ f + °°j jAt ^ 18M7r3 

A i = TFT / dz hi*) 



NA ' 35AN ' 

vr 3 f + °° „. . „. . 2Mtt 3 

Herein the ultralow-energy fields 4>(x) and /i(x) only have been retained, whilst the heavy 
scalars and fermions with masses ~ M have been decoupled. 



3 To be precise, the ultralow-energy effective Lagrange density for the light states on the four- dimensional 
Minkowski space-time is well defined up to subtraction of an infrared divergent constant. 
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Quite remarkably, the domain wall Lagrange density ( 4.23 ) has a non-trivial large cut- 
off limit provided that the ratio M/A < 1 is fixed. The four-dimensional ultralow-energy 
theory happens to be interacting with the ratios, 



gf : Ai : A 2 : A 3 ~ 6 : 5 : 4 : 3 , (4.24) 

being independent of high-energy scales and regularization profiles - here we leave aside 
the issues concerning the renormalization group improvement. 

However the solution (J) is not of our main interest because the vacuum expectation 
value of the field H vanishes and does not supply the domain wall fermion with a light 
mass. 

5. Domain walls: Higgs phase 

Evidently, the domain wall solution of Eq. (|4.1| ) as well as the constant background solution 
of Eq. ( |3.7D just break the r\- and T2-symmetries of the Lagrange density ( |2.2| ), whereas 
they keep the T3-invariance untouched - see Eq. ( |2.5| ) . Meanwhile, the second domain wall 
background of Eq. ( |4.2[ ) does break all the r-symmetries, i.e. it realizes a different phase in 
which the masses for light particles are naturally created. We notice however that for the 
kink solution {K) realizing the space defect in the z-direction the combined parity under 
the transformations, 

z — ► — z, &(x,—z) — ► — <&(x,z), H(x,—z) — >H(x,z), 

$i(x, -z) — ► 7 5 r 3 ^/(x, z), ^i(x, -z) — > ~^>i(x, z)%t 3 , (5.1) 

remains unbroken. 

As we will see below the mass scale for light particles is controlled by the parameter 
H = V2A2 — M 2 , which describes the deviation from the critical scaling point where the 
two regimes of Eq.s ( J4. 1| ) and Q4.2j ) melt together. 

As we want to supply fermions with masses much lower than the threshold of penetra- 
tion into the fifth dimension, to protect the four-dimensional physics, we assume further on 
that \x <C M. At the solution (K): S(X) = Sk(z), the 2x2 matrix-valued mass operator 
M[S'(X)] of Eq. ( |4.6| ) for scalar excitations gets the following entries, 

M^ 1 = -d 2 + 4M 2 + 2{fi 2 - 3M 2 ) sech 2 ((3z) , (5.2) 
Mf = -d 2 + M 2 -[i 2 - 2(M 2 - 3fi 2 ) sech 2 (/3z) , (5.3) 
= Mf = AMn sinh(/?z) sech 2 (/?z) , (5.4) 

the positive quantities fi and being defined in Eq. (|4.4| ). As this mass operator is non- 
diagonal it mixes the scalar fields $ and H. However, this mixing does fulfill the combined 
symmetry, 

s(x, z)=±a 3 s(x, -z) , M[S K (z)} = a 3 M[S K (-z)] a 3 , (5.5) 

because the diagonal elements of M[S#-(z)] are even and the off-diagonal ones are odd with 
respect to the reflection z — > —z. This symmetry allows to classify the non-degenerate 
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eigenstates of the operator Mx(^) = M [Sk(z)] according to their parity. Actually, they 
may consist of even-odd or odd-even pairs of components 



S(X) - S K 



si(x,z) 

k S2(X,Z) t 



One more symmetry can be revealed in respect to the simultaneous change in the sign of 
\i and the reflection z — > — z, that means 



M[S K (z,fi)]=M{S K (-z,-n)) . 



(5.6) 



This formal symmetry turns out to be crucial in perturbation theory to build up eigenstates 
of ~M.k(z) for small [/, <C M. Indeed, after a proper normalization of the eigenstates, the 
expansion in powers of [i of the latter ones is constrained as follows: the parity-even 
components of s(x, z) have to contain even powers of /i, whilst the parity-odd ones have to 
contain odd powers of fi. 

As a matter of fact, for fi <C M the operator Mx(z) can be suitably decomposed 
into its diagonal part with the very same diagonal matrix elements of Eq.s ([O]) and 



(4.10), in which M is replaced by (5 and [i = 0, plus the small perturbation AM/c, 



M K = + AM A - 



(5.7) 



where 



whereas 



M 



(0) 
K 



-d 2 + 4/? 2 - 6/3 2 sech 2 (/?z) 




-d 2 z + (3 2 - 2/? 2 sech 2 (/3z) ^ 



AM 



K 



4(3 2 e 



etanh 2 (/?z) 
v 7 ! + e 2 tanh(/3z) sech(/3z) 



Vl + e 2 tanh(/3z) sech(/3z) 
e sech 2 (/3z) 



Herein the dimensionless parameter e = n/ f3 <ti 1 controls the deviation from the scaling 
point where the two regimes (J) and {K) do coincide and both scalars are massless. One 
can choose the parameters (3 and e as independent ones, characterizing the mass unit 
and a scaling-point deviation. Then the initial mass scales M and \i can be conveniently 
expressed in terms of the latter ones as 



M = (3^fi 



+ e 2 
r(0) 



H = /3e . 



(5.8) 



Now, it turns out that the operator has two zero-modes, as it immediately follows 
from the comparison with the related operators (4.9) and ( [4.10 ) in which M — > /3 and 
M 2 — 2A2 is omitted. The corresponding massless eigenfunctions are obviously given by 
Eq.s ( 4.15| ) and ( 4.16| ) up to the replacement M — ► /3. The deviation from the scaling point 
is therefore fully generated by the perturbation term AMx, as it does. 

The masses of the lightest localized scalar states can be obtained again - see Eq. (|4ll| ) 
- from the Schrodinger equation with a matrix-valued potential: namely, 



£ M%M K (X) = m 2 W K (X) 



1,2. 



(5.9) 



=1,2 
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It is easy to find the massless state as a zero- mode of the operator in Eq. (|5.7|) . Indeed, 
differentiating with respect to z the equations ( [4.3| ) of the stationary configuration, one 
obtains the normalized zero-mass solution from the derivative of the kink-like solution 
Sk(z), with the upper signs, in the form 



V 

7x 



3/3 
4 + 6e 2 



= d z S K (z 



/? 3 (4 + 6e 2 

v / TTe 2 " sech 2 (/3z) 
e sinh(/3z) sech 2 (/32:) 



(5.10) 



according to the translation symmetry breaking in the kink background [47]. Evidently, 
this eigenfunction is related to the unperturbed zero-mode of Eq. (|4.15| ) at \i = which is 
parity-even. Respectively, it consists in turn of the parity-even upper component 4^{z) 
and the parity-odd lower component (j)J,(z), in accordance to the combined symmetry 
of Eq. ( |5.5|) , because the perturbation term does not change the parity properties of an 
eigenstate. Finally, one can convince oneself that, for a given /3, the expansion in e *-* fx 
contains only even powers of this parameter for the upper, even component and only its 
odd powers for the lower, odd component in accordance with the /x-symmetry of Eq. (|5.6|). 

The second light scalar state O^-(X) arises from the perturbation of the zero-mass 
state of Eq. ( 4.16| ) for fi = 0. Therefore, it contains a parity-odd upper component h~^{z) 
that can be expanded into odd powers of e, together with a parity-even lower component 
hX,(z) which can be expanded in even powers of e, in such a way that the exact normalized 
second light scalar state can be written as 



n h K (x) = h(x) 



h K (z) 



A 



(5.11) 



To the first order in the e expansion - see Appendix C - the localization functions read 



h K {z) 



h K {z) 




( h + K {*) . 


V 



-e(3z sech 2 {(5z) 
sech.{f3z) 



+ G(e 2 



The mass of this scalar state, to the first order in e expansion, is given by 

m h = (3eV2 + 0{e 2 ) ~ . 



(5.12) 



(5.13) 



One can see that now the mass eigenstates enter both in the upper component si(X) = 
<3?(AT) — (<3?(X))o and in the lower component of S2(X) = H{X) — (H(X))q of the low-energy 
scalar field: namely, 



S1 (X) 

S2(X) 



(x)<f>+(z) + h(x) h K (z) 
(x)cf> K (z) + h(x) h+(z) 



A 



A 



(5.14) 
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However, the admixture of the opposite-parity states is strongly suppressed in the ultralow- 
energy effective Lagrange density in the Minkowski space-time, owing to the integration 
over the extra-dimension and to the high order in e. 

As well as for the solution (J), in the vicinity of the vacuum solution (K) of Eq. (|4,2|), 
the ultralow-energy effective Lagrange density for the light states can be obtained from 
Eq. ([O]) in a similar form. In the leading and next-to-leading order of e-expansion, one 
can show that the only difference concerns the appearance of the mass for light fermions and 
of a cubic scalar interaction that we calculate for positive values of {&(X))q and (H(X))q 



in Eq. (L2): namely, 

6$ [?,^,<M] =4 4) -m/?(x)V(x)-A 4 /i 3 (x)+0(e 2 ) , (5.15) 

the fermion mass being determined by the expression 

f +OD - TT 7T 

771/=/ dzip (z)H K (z)i/}o(z) = -fie = -/i , (5.16) 



where Eq.s (|4.2| ) and ( 4.17 ) have been employed, whereas 



8tt 3 
AiV 



-oo 



dz 

oo 



2 



{^ K (z) h- K {z)+H K {z)h + K {z)} 



Mir 

"Han- (5 ' 17) 

We stress that a generally a priori possible 3-point vertex 4> 2 {x) h(x) does actually receive 
mutually compensating contributions from different mixing terms. Thereby the direct 



decay of the massive Higgs-like boson h into a pair of massless branons [24] is suppressed 
and the low-energy Standard Model matter turns out to be stable. 

To sum up, in the present model the ratio of the dynamical fermion mass (presumably 
the top quark one E3j) to the Higgs-like scalar mass is firmly predicted to be 



m h : nif 



4>/2 : 7T ~ 1.8 , (5.18) 



which is somewhat less than such a ratio as predicted in the Top-mode Standard Model - 
the Nambu relation gives rrih : mt = 2. 

Finally, concerning the coupling parameters g 2 , X\, A2, A3, they are essentially de- 
scribed by Eq.s (|4.23| ), up to the two orders in e-expansion, the corrections being of the 
order 0{e 2 ). We would like to stress that the coupling of fermions to the massless scalar 
4>(x) does no longer appear even in the next order ~ e 2 , because its mixing form factor 
4>x(z) in Eq. ( |5.14 ) is an odd function, making thereby the relevant integral of Eq. ( 4.23| ) 



identically vanishing. 



6. Manifest breaking of translational invariance 

One can conceive that in reality the translational invariance in five dimensions is broken not 
only spontaneously but also manifestly due to the presence of a gravitational background, 
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of other branes etc. In a full analogy with the Goldstone boson physics one can expect [24] 
that the manifest breaking of translational symmetry is such to supply the branons with a 
mass. 

In the model presented in our paper the natural realization of the translational symme- 
try breaking can be implemented by the inhomogeneous scalar background fields coupling to 
the lowest-dimensional fermion currents. Let us restrict ourselves to the scenario of the type 
(K) and introduce two scalar field defects with the help of the functions F$(z) and Fh(z), 
which are supposed to be quite small, i.e., \F$(z)\ <C |(<3?(X))o| and |F#(z)| <C \(H(X))q\. 
These scalar defect fields catalyze the translational symmetry breaking and the domain 
wall formation by means of their interactions with the fermion currents, 

4 5) = -F$(z)^(X)r 3 *(X) - F h (z)^(X)t^(X) . (6.1) 



When supplementing the Lagrangian ( |2.2j ) under the Hubbard-Stratonovich transforma- 
tion, one can see that those background defect fields do actually couple to the auxiliary 
scalar fields in Eq. ( |2.2D after the replacements 

$ , — > $-F 4 (z) , H i — ► H-Fh(z) . (6.2) 

The latter redefinition of the auxiliary fields reveals the explicit coupling of defect fields to 
the auxiliary scalars, which dictates in turn the eventual change of the scalar field Lagrange 



density (3J): namely. 



£(5) [q,y,<!>,Fz(z),H,F H (z)] = £® + 

2NA 3 . . 2NX 3 , . NX 3 ^ 2 , NX 3 , , , , 

$F$(z) + HF H {z) ¥%{z) F 2 H {z) . (6.3) 

9i 92 9i 92 

In order to prepare the scalar part of the low-energy effective action in comparable units 
let us re-scale, 

- £^2 W . F h(z) ee ^ Hz) , (6.4) 

where we have kept in mind that 9i — 92 — 9 CT = 9tt 5 . If we assume the dimensionless 
functions f$ and f# to be O(l), then the last two terms in Eq. (|6.3[) are of the order fi/A 
and thereby negligible within the approximation adopted in this paper. 

Taking the notations of Eq. ( [4.4D into account, it is not difficult to show that the 
effective Lagrange density for low-energy scalar fields in the presence of defects can be cast 
in the form 



+ (d Q $) 2 + (d a H) 2 + 2Af 2 $ 2 + (M 2 + fi 2 )H 2 - ($ 2 + H 2 ) 2 ] . (6.5) 



The stationary vacuum configurations of scalar fields obey now the modified set of Eq. (^) : 
namely, 

d a d a <5> = v 3 h(z) + 2$ [ M 2 - $ 2 - H 2 ] , 

d a d a H = ii 3 i H {z) +H [ M 2 + fi 2 - 2H 2 - 2$ 2 ] . (6.6) 
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Let us search for the solutions of these equations generating a domain wall of type (K) 
which preserves the symmetry (|5.1|) . Evidently, the latter ones can be actually achieved if 
both the scalar defect fields and the very vacuum solutions have the same definite parity 
properties, i.e., i$(z) and <&k(z) being odd functions of z, while £h(z) and Hk{z) being 
even ones. Thus we tune ourselves to what we could call the collinear mechanism of a brane 
creation. As we shall see here below, it allows to formulate within the perturbation theory 
the self-consistent catalyzation of a domain wall in the presence of some weak defects, i.e., 
for sufficiently small defect functions f$(z) and {h(z). 

Indeed, let us accept the latter assumption and write down the trial solutions as the 
sum of the unperturbed functions &k(z) = Mtanh(/3z), Hk{z) = fJsech(Pz) given in 
Eq.s ( [4.2D and the corresponding small deviations 5f$>(z),5fH(z) due to the presence of a 
weak defect, so that 



fU [tf ( f*) J (h k (z) + 6 { H(z] 
To the first order in deviation functions, Eq.s fl6.6|) entail 

M[S K (z)} 5 f S = M K (z)6 { S = fiH , f{z) 



S K (z) + 6 f S(z) . 



h(z) 
{ 1h(z) 



(6.7) 



(6.8) 



where the second variation operator M^(z), already introduced in Eq. (|5.7|) , does appear 
in this inhomogeneous equation. Its solution can be suitably presented by means of the 
spectral decomposition for the operator ~M.k(z) itself. In fact, as it was elucidated in the 
preceding Section, its spectral resolution contains two sets of eigenfunctions with opposite 
parity properties. In particular, once we search for solutions of odd upper component 
and even lower component type, the relevant part of the spectral decomposition of the 
operator M^(z) consists of the projector onto the bound state, Eq. ( |5.12| ), belonging to 
the eigenvalue m\ ~ 2// 2 , together with the projectors onto the continuum improper states, 
starting from the eigenvalue M 2 3> 2n. 2 : namely, 



M K (z) 



m 2 h P h + 



M 2 



da g(a)a P CT + irrelevant 



(6.9) 



where P/,,P ff are the projectors onto the proper and improper states with the relevant 
parity properties respectively, whereas a denotes the invariant mass of the continuous part 
of the spectrum and g(a) the corresponding density of the states. After inversion of the 
operator (^) the approximate solution of Eq. fl6,8|) reads 

,2 



5fS ~ [i 



da g(a) — P CT f 

M 2 G 



(6.10) 



Thus, taking into account that u 2 <C M 2 < a, a finite solution ~ \i exists iff the projection 
of the defect function f(z) onto the bound state is not vanishingly small. Then the first 
order approximate solution is given by 



StS(z) 



^^h K (z) + 0{fie 2 ,^ 2 ) , 



(6.11) 
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where the dimensionless parameter £ is introduced by means of the relation, 



dz 



h K (z)h(z)+h + K (z)i H (z) 



(6.12) 



Let us first assume that e <C |£| <C 1. Then, for a wide variety of scalar defects f(z), the 
dynamical mechanism of the domain wall formation is essentially triggered by the light 
bound state of the mass operator M#, up to the first order approximation in perturbation 
theory. Hence, for this variety, the manifest breaking of translational invariance is labeled 
by a one-parameter family of weak defect amplitudes |£| < 1. However, it turns out that 
in the calculation of the branon mass the next-to-leading terms - quadratic in £ - are indeed 
required. The latter ones are more sensitive to the high-frequency modes in the spectrum 
of the operator Mjj-, owing to the non-linear interaction in Eq.s (|6.6| ). Consequently, they 
appear to be dependent upon the details of the defect function i(z). The analysis concerning 
the general solution will not be addressed here, instead of, let us make the ansatz for f(z) 
which allows to solve Eq.s ( |6.6[ ) exactly. Firstly, one is evidently free to restrict the allowed 
set of defect functions in such a way that the solution ( |6.11 ) keeps itself reliable all over 
the range e <C |£| < 1. In particular, it is provided by the ansatz (A) : 



6^^(<5 f ^) 2 + 4H K 5{H 5 { <f> + 2{5 { $) 3 + 2(5 f #) 2 <5 f $ ; 



a¥ A) 



m K {5fHf + <5 f $ 5 t H + 2(5 { Hf + 2(5 f <S>) 2 5 i H , 



(6.13) 



where Hx are given by Eq.s ( |4.2| ) and the deviations 5f$, 5{H are taken from Eq. ( |6.11| ). 
For such a defect, to all orders in the weak defect amplitude, the modified kink configuration 
of the scalar fields reads 



! - 2 ^ tanh(^) - i< 2 
e(l + f ) sech((3z) 



[1 + £e*) tanh(/3z) - £e 2 /3z sech 2 (/3z) 



(6.14) 



As the lightest domain wall fermion gets its mass only from the lower component of S^(z), 
then the mass correction is just reproduced by the following multiplicative factor of the 



unperturbed mass (|5.16[) 



7T 



m 



f 



(6.15) 



In terms of j3 and e, the mass operator for scalar particles receives an additional perturba- 
tion as compared to Eq.s ([5/1]): namely, 



M 



(-1) 
K 



M.P + AM K + AM A 



(6.16) 



AM A = 4/3 e£ sech 

e + ~3(3z tanh(/3z) 



x 



X 



sinh(/3z) 



s'mh((3z) 
e [3+*£-p z tanh(/3z)] ^ 



[l + 0(e 2 ) 
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where the basic definition of the second variation operator Eq. ( |4.q ) has been used. There- 
from one obtains the masses of scalar states - see Appendix D. In particular, the lightest 
branon is no longer massless 

(6.17) 



and its localization function turns out to be 



ye 

5 



<pa{z) ~ y t 



sech 2 (/3z) 
k -e (1+0 sinh(/?z)sech 2 (/?z) ^ 



Respectively, the correction to the Higgs boson mass reads, 



4^)W(i + 3e + ^ 



whilst the corresponding localization function is given by 



h A (z) 



-e(l + f) l3zsech 2 (/3z) 
sech(Pz) 



(6.18) 



(6.19) 



(6.20) 



We see that the regular defect does not rigidly polarize the domain wall vacuum (K) in 
a fixed direction - the sign of £ can be not only positive but also negative while delivering 
a local minimum of the induced scalar field Hamiltonian, i.e. physical masses for scalar 
bosons. For the latter case, £ < 0, one can reduce the mass of Higgs particle as compared 
to the fermion (" top-quark") mass, 



ni 



(A) 



m} 



32(1 + 3g + 5£ 2 /3) 

7T 2 (l+0 2 



(6.21) 



For instance, if £ — —0.4 and the fermion mass is assumed to be of order of the top-quark 
mass ~ 175 GeV, then the Higgs mass is estimated to be ~ 135 GeV, which is acceptable 
from the phenomenological viewpoint [ja2]j . For the same choice the branon mass is found 
to be ~ 100 GeV. However, we stress that the predictions for scalar masses are essentially 
based on contributions which are quadratic in £, especially for the branon mass. But the 
latter terms depend on a model for the space defect. 

Both the localization functions (3.18), ( 6.20| ) contain corrections linear in £ which are 
the only ones relevant for calculation of the induced 3-point and 4-point vertices for low- 
energy scalar fields. The coupling constants g 2 ,Ai,A2,A3 in the 4-point vertices on the 
domain wall are built of the components (f)\,h\, which are universally given by Eq.s fljjj ) 
at the leading order in e 2 and thereby do not depend on a small background defect. The 
coupling constant A4 in the 3-point vertex h 3 (x) is described by the integral in Eq. ( 5.17| ) 
which is evidently homogeneous in the multiplier (l + £). Hence, in the presence of a defect, 
this coupling is simply renormalized as follows, 



:i+oa 4 ~ (1+omt 



Mvr 

In 



(6.22) 
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Another possible induced coupling constant in the vertex (j) 2 (x)h(x) is also homogeneous 
in the multiplier (1 + £). Thus it remains vanishing, the Higgs boson does not decay into 
branons and the Standard Model matter is stable. 

The previous ansatz (A) has been a representative of non-top ological or regular defects 
with background functions decreasing at the infinity. Now we shall introduce and discuss 
another kind of defect of a topological type, the influence of which on the branon masses 
is drastically different. Let us take the following ansatz (B) for the defect i(z) providing 
the analytical solution for vacuum configurations, 



^ B \z 



(j, Ksech(/3z) 



(6.23) 



where 7, n are dimensionless parameters. Notice that the components are just borrowed 
from the unperturbed solution (JO) . We accept that the upper component can be enhanced 
including the normalization on M rather than on fi - otherwise it becomes irrelevant as it 
will be clear later on. We also remark that, on the one hand, the related function is not 
square integrable and lies outside the domain spanned by the normalizable eigenfunctions 
of the operator M^. On the other hand, it represents a genuine topological defect whereas 
the lower, square integrable component does not influence the field topology, i.e., it is a 
regular defect. 

The exact solution of Eq.s ( |6.6D keeps the form of Eq. Q4.2D , although its coefficients 
are now controlled by 7, k: namely, 



M atanh(/3z) 
H (1 + £) sech(^z) 



(6.24) 



The constants a, £, (3 obey the relations, 

/2 



2(a 2 



Da 



7£ 



1 + e 2 



1 + 



2£ + £ 2 + ^) (l + = «, 



M 2 a 2 - fi 2 (l + 2 ~ P 2 



l + € 2 [t 1 -2i-e 



(6.25) 
(6.26) 

(6.27) 



From the last equation one can reveal a certain interplay between the enhanced topological 
defect and the regular defect. Instead of the inversion of Eq. ( 6.26 ) we prefer to use £ as 
an independent parameter, thereby applying this relation as a definition of n. 
Let us expand now the solution (|6.24| ) in powers of e, 



S { k\z) 



(6.28) 







1 + |e 2 (2 + 7) 1 tanh(/3z) + ±e 2 ( 7 - 4£ - 2£ 2 )/3z sech 2 (/3z) 
e(l + £) sech(Pz) 



Here one can see that the linear analysis in Eq. ( |6.11[) is again valid and entails £ 



(k/2) — (7/4) in accordance with the linear part of the exact solution (6.26). As the lower, 
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even component of Eq. ( 6,2§| ) coincides with that one in Eq. ( 6.14| ) , the shift in the fermion 



mass remains the same as in Eq. ( 6.15Q 



The mass operator for this ansatz (B) gets a similar functional structure albeit different 
coefficients: namely, 

) = + AM K + AM B , 
AMg 1 = j3 2 e 2 sech 2 (/3z) [4£ + 2£ 2 

+3(7 -4f -2^ 2 )/3ztanh(/?z) +37sinh 2 (/3z)] + O (e 4 ) , 
AM 1 ? 2 = AMf = 4£/3 2 esech(/3z)tanh(/?z) + O (e 3 ) , 
AMf = /? 2 e 2 sech 2 (/3z) [l2£ + 6£ 2 

+(7-4^ -2e 2 )/?ztanh(/3z) +7sinh 2 (/?z)] + O (e 4 ) . (6.29) 

The corresponding mass spectrum of lightest scalar states is calculated in Appendix D. In 
particular, the branon mass is no longer governed by terms quadratic in £, 

mf ) ) 2 «7M 2 - (6-30) 

Now, the striking evidence occurs for a strong polarization effect induced by a topologi- 
cal defect: the local minimum is guaranteed only for asymptotics at infinities which are 
coherent in their signs, i.e., for positive 7. 

The Higgs mass encodes both topological and non-top ological vacuum perturbations, 



m i B) ) 2 ~ m 2 ( 2 + 6 £ + 3 ^ 2 + \y) ■ ( 6 - 31 ) 



As 7 > the topological defect makes the Higgs particle heavier. However, once again, the 
sign of £ is not fixed by the requirement to provide a local minimum. Therefore one can 
find a window for a relatively light Higgs scalar. If branons are very light, viz. £ 3> 7, then 
one can neglect the last term in Eq. (3.31) and reproduce the Higgs mass ( |6.19| ) of ansatz 



(A) with a reasonable precision. 

As to the localization functions for scalar states, they coincide with those ones of ansatz 
(A). Thus they are universally parameterized by terms linear in £, to the leading orders 
in e-expansion. Respectively, the induced coupling constants g 2 -, Ai, A2, A3, A4 appear to be 
universal as well. 



7. Conclusions and discussion 

The main task of the present work has been the explanation of how light fermions and 
Higgs bosons may be located on the four-dimensional domain wall representing our Uni- 
verse. The five-dimensional fermions with their strong self-interactions have been used to 
discover the vacuum that breaks spontaneously the translation invariance at low energies. 
Its remarkable feature is the appearance of a topological defect which induces the light 
particle trapping. Depending on the relation of the four-fermion coupling constants, the 
light fermion states may remain massless or achieve the masses ~ fi, supposedly, much less 
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than the characteristic scale M of the dynamical symmetry breaking. These two phases 
differ in the level of symmetry breaking. Namely, the internal discrete r-symmetry of the 
initial fermion Lagrange density is generated by the algebra of Pauli matrices and this sym- 
metry can be broken partially or completely. In the latter case, the light fermions acquire 
a dynamical mass which is assumed to be ~ fi <C M and the ratio of the Higgs boson mass 
to this fermion mass is ~ 1.8 close to the ratio in the Top-mode Standard Model p0| ]. 

As a consequence of spontaneous breaking of the r-symmetry and of translational 
invariance two composite scalar bosons arise, one of which is certainly a Higgs-like massive 
state whilst the other one is a Goldstone-like state - a massless domain-wall excitation 




called branon [23, 24]. On the one hand, the latter particles turn out to be sterile because 
in the ultra- low energy world they do not couple directly to the light fermions. On the other 
hand, their original coupling to fermions in Eq. ( |3.lD mixes light fermions and the ultra- 
heavy fermion states with masses ~ M. Evidently, the amplitude of branon production in 
the light fermion annihilation due to exchange by an ultra- heavy fermion state ~ g 2 m 2 /M 2 
(see Fig. 1) is strongly suppressed if the domain wall thickness ~ 1/M is much less than 
the Compton length for light particles, 1/M <C 1/rrif. 

However, it is quite plausible that the f ^ 
translation and r-symmetries are manifestly 
broken just because the five-dimensional world 
contains a classical background of gravitational 
and various matter fields with a net result 
similar to a crystal defect. In the present 
paper we have assumed the scalar nature of 
a five-dimensional space-time defect, which 
triggers the domain-wall formation through 
its coupling to the low-dimensional fermion 
operators. The two types of such defects have 
been investigated. 

The first type is a non-top ological one, the shape of which belongs to the Hilbert space 
of square integrable functions. In this case it happens that the only relevant dimensionless 
parameter £ governing the manifest symmetry breaking is given by a projection of the 
defect profile function on the Higgs boson wave function. The occurrence of such a defect 
supplies the branons with a small mass ~ £/i though its numerical value is rather model 
dependent. Meantime, the polarization of the domain wall background, i.e. its sign, is not 
correlated with the sign of a non-topological defect - for any sign, the local minimum of 
the low-energy effective action is attained. 

The second type of defects is a topological one, with different non-vanishing asymp- 
totics at the infinity. It makes much stronger the catalyzation of a domain wall, namely, 
the minimum and respectively the physical branon masses are reached only for coherent 
signs of the profiles of defect and domain-wall. In the latter case, again, the branons are 
supplied with a small mass. 

Concerning the Higgs boson masses, their ratio to the fermion (~ top-quark) mass 
can be substantially reduced with an appropriate choice of a non-topological part of the 



f Y 

Figure 1: Branon (j> creation from light 
fermions /, / due to super-heavy fermion ^ ex- 
change. 
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defect ~ £. In particular, the Higgs masses may be well adjusted to a phenomenologically 
acceptable [|5^] value ~ 135 GeV for a reasonably small value of a defect £ ~ 0.4. 

In all scenarios the ratios of the induced coupling constants in the four-dimensional 
effective action for light particles are predicted unambiguously and in a model-independent 
way. Their absolute values are universally parameterized by the ratio of the symmetry 
breaking scale M to the high-energy (compositeness) scale A - a cutoff in the four-fermion 
interaction. This ratio M/A < 1 can be thought of as an empirical input which must be 
established from the relevant experiments by measuring, for instance, the Yukawa coupling 
constant gj in the direct Higgs production via the fusion mechanism, 
j ± However, the branon matter is much more elu- 

sive for experiments because it is not produced by 
fermions at the leading order in the Yukawa cou- 
pling and does not appear as a decay product of 
the Higgs bosons. Still the fermion annihilation 
(see, Fig. 2) may result in missing energy events 
due to radiative creation of branon pair mediated 



f 9 by Higgs bosons. The corresponding Feynman in- 

tegral is well convergent and therefore, at very low 
Figure 2: Missing energy in fermion an- momenta; the amp litude of this process can be esti- 

nihilation due to branon emission. , , , , r , o\ 2/2 m • i.i 

mated to be of order gj ■ M- m^/fi. We notice that 

if they are produced by very light fermions, e.g. e + e~ , the very last ratio is negligibly small 

putting such signals into the experimental background. Nonetheless, for top-quarks (again 

in the fusion production) there might be a possibility to discover branon pair signals for 

sufficiently light branons. More qualitative estimates will be presented elsewhere. Anyhow, 

the sterile branons seem to be good candidates for saturation of the dark matter of our 

Universe - see Ref. [EJl for a more detailed motivation. 
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A. One-loop effective action 

Let us consider the second order matrix-valued positive elliptic differential operator 

D*D = -d„d„ - dl + M 2 {X) = -d 2 + M\X) , 
M 2 (X) = $ 2 (A) + H 2 (X) - r 3 fl$(X) - n ftH(X) , (A.l) 
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acting on a ra-dimensional flat Euclidean space (eventually n = 5). Our aim is to evaluate 
the matrix element of the distribution of the states operator: namely, 

(X\$(Q 2 ~ tfD)\Y) = [ C+l °° ^ eXp{ l Q2} (X\ exp{-tD^D}\Y) , c> . (A.2) 

As a matter of fact, if the positive operator D'D is also of the trace class, then the trace of 
•d{Q 2 — D^D) does represent the number of the eigenstates of D^D up to the momentum 
square Q 2 . It is convenient to write the heat kernel in the form 

(X\exp{-tD^D}\Y) = (47rf)-"/ 2 exp j - (X ~ Y) * j Sl(t\X, Y) , (A.3) 

where £l(t\X, Y) is the so called transport function which fulfills 

d t + ^-+ri x D x )si(t\X,Y) = 0, (A.4) 
]hnQ(t\X,Y) = 1 , (A.5) 

in such a way that 

lim (X\ exp{-tD^ D}\Y) = 15^{X - Y) . (A.6) 
If we insert Eq. ( |A.3 ) in Eq. ( A.2| ) and change the integration variable we obtain 

{X\d{Q 2 - D^D)\Y) = 2Q n (Air)- 1 - n / 2 / dt r 1 ^ 2 

J c—ioo 

Now, if we write 

/ t \ M / t \ 

n^\X,Y) = Y)Q- 2k + n [n/2]+l \X,Yj , (A.8) 



it turns out that, by insertion of Eq. ( A, 8 ) into the integral ( |A.7| ), the last term in the RHS 
of the above expression becomes sub-leading and negligible in the limit of very large Q. As 
a consequence, the leading asymptotic behavior of the matrix element ( |A.2| ) in the large Q 
limit reads [53] 

[n/2] 

(X\tf(Q 2 - D*D)\Y) (4vr)-"/ 2 ^ a k {X,Y)Q n ~ 2k 

k=0 

c+lc °± t ^-n/2 L Q 2 (x-y) 2 

c-ioo 2m \ At 

v v [»/2] n n-2k 

fc=0 v ' ' 
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In the diagonal limit X = Y, for n = 4, 5 the relevant coefficients of the heat kernel 
asymptotic expansion take the form 54] 



a (X,X) = l; ai (X,X) = -M 2 (X) ; 

a 2 (X,X) = \[M 2 {X)] 2 - -d 2 M 2 {X) , (A.10) 
2 o 

in such a way that we can finally write the dominant diagonal matrix element - leading 
eventually to the five-dimensional Euclidean effective Lagrange density - as 



{X\$(Q 2 -D^D)\X) 

Q- 2 M 2 (X) + ^Q- 4 [M' 2 (X)] 2 - ^Q- 4 d 2 M 2 (X)} . (A.ll) 



60vr 3 \ 2 

B. Spectral resolution for Schrodinger operators 

In this Appendix we calculate the spectra and eigenfunctions of the Schrodinger operators 
Mj 1 and M 22 arising in Eq.s (|4.9| ) and ( [4,10 ) and presented in the factorized form in Eq.s 



( |4.12[) and (|4.13|) . Let us choose the inverse mass units, z = y/M and introduce the 



dimensionless operators 

M 1 / = M 2 4 = M 2 q+q- , (B.l) 

Mj 2 - M 2 + 2A 2 = M 2 d\j = M 2 q+ qjj , (B.2) 

9* = Td y + 2tanh(y) , (B.3) 

Qh = Td y + tanh(y) . (B.4) 

The two operators and can be embedded into the super symmetric ladder when one 
adds their supersymmetric partners 

= <h, <lt = d% + 3 , d H = q H qj I = -d 2 + 1 , (B.5) 

the latter one describing a free particle propagation. The intertwining Darboux relations 
read 

5$ d% = q^ , 
q± dfj = d%q± , 

q~H 9$ <4 = qjjidjj + 3)q$ = (-d 2 + A)q]j q^ , 
d% q% Qh = 9$ + 4 ) - 

iH d H = (~ d y + !)?H ' 

d + H q + H = q + H {-d 2 y + l) . (B.6) 

From these relations one concludes that the spectra of the operators dt and djj are almost 
equivalent one to each other and, in turn, both equivalent to the free particle continuous 
spectrum - up to a shift of the origin. In other words, they consist of the continuous part 
and of the zero-modes of the intertwining operators q~j^q^ and qjj respectively. Owing to 
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Eq.s 
namely, 



their eigenvalues and normalized improper eigenfunctions are related as follows: 



ipq>(ky) 

ipH{ky) 
i>o{ky) -- 



^k 2 + 4 
QH^oiky) 



Vk^Ti 



v / (fe 2 + l)(A; 2 + 4) 

Em 



E%(k) = k 2 + A 



k 2 + 1 



1 



exp(iky) 



k£R 



(B.7) 



The imprecise equality = signifies a possible discrepancy due to the existence of zero- modes. 
In particular, for the operator d~jj the continuum starts from -E^(O) = 1 and the only zero- 
mode ho(y) = (l/\/2) sech(y) just coincides with the zero-mode of qjj. Concerning the 
operator dt, one can easily verify that the continuum starts from E^(0) = 4, whereas the 
zero-mode 4>o(y) = (\/3/2) sech 2 (y) coincides with the zero-mode of gj. Furthermore, the 
second parity-odd proper eigenstate 



q$ h o(y) = vW 2 " tanh(y) sech(y) 



(B.8) 



does coincide with the zero-mode of q H q^ and corresponds to the first non- vanishing, 
discrete eigenvalue E$ = 3. 

C. Perturbation theory for the first excited state 

In this Appendix we develop the perturbation theory for the first excited eigenstate and the 



corresponding eigenvalue of the operator M# defined in Eq. (5.7). As in the Appendix B 
the inverse mass units z = y/(3 are chosen, so that the dimensionless operator 3 = M^//? 2 



has the following components in terms of the operators in Eq.s (B.1)-( B.4 ): namely, 

e sinh 2 (y) 



f4 







d 



H 



+ 4e seclr(y) 



(C.l) 



V / T+~e 2 sinh(y) 
VI + e 2 sinh(y) e 

The second, perturbation term contains the elements of different order in e. This is the 
ultimate reason why, in order to compute the corrections to the eigenvalues, we have to take 
into account not only the eigenfunction of zeroth-order, but also the first-order correction 
turns out to be necessary. Thus we use the eigenfunction hx{y) = [hj^(y), h^ < (y)] from 
Eq. ( |5.12 ), taking into account that, to the lowest order, the parity-odd function h~^-(y) is 
of order e and the parity-even one is h^(y) ~ ho(y) = (1/^/2) sech(y). To the first order 
in e, the function h~^(y) is derived from the following equation 

(y) = — 4etanh(y) sech(y) /io(y) = -= tanh(y) sech 2 (y) 



V2 



-j=q^sech 2 (y) , 



which follows from the basic definition E/ir- = E^hx and the estimation of E^ ~ e 2 
there are no normalizable zero-modes for the operator q~^ one concludes that 



h K (y) 



-^sech 2 (y) 



(C.2) 
Since 

(C.3) 



- 28 - 



which represents a first-order differential equation with the normalizable solution 



h x(y) = --^2/sech 2 (y) 



(C.4) 



Now, the eigenvalue Eh can be obtained from its integral representation Eh = (hx | S | hx) 
in terms of the normalizable function (h^hx) = 1- To the first order in e it reads eventually 

n + OO 

L), 1 / dy ho(y) (e 2 sech 2 (y) h (y) + etanh(y) sech(y) h K (y)} = 2e 2 . (C.5) 



D. Perturbation theory in the presence of defects 

The inclusion of a small regular defect is described by the scalar functions f$(z), £h( z ) and 
leads to the change of the mass matrix M#- presented in Eq. (6.16) for the ansatz (A) and in 
for the ansatz (B). If we introduce the dimensionless operators Sj , j = A,B , 

we obtain 



Eq. (ET2 

then according to the notations of Eq. (| 



Ej = E K + Aj , 
A 11 



A 12 

A 22 
^A 

A 11 

a 12 

a 22 



e 2 sech 2 (y) [4£ + 2£ 2 - 12y tanh(y) ] , 
A 2 / w 4e£ sech 2 (y) sinh(y) , 
e 2 sech 2 (y) [ 12f + 6f - Ay tanh(y) ] , 

e 2 sech 2 (y) [ 4£ + 2f + 3(7 - 4£ - 2£ 2 ) y tanh(y) + 37sinh 2 (y) 

A 21 
^B 



4^e sech 2 (y)sinh(y) , 
e 2 sech 2 (y) [ 12f + 6£ 2 + (7 - 4£ 



(D.l) 
(D.2) 

(D.3) 



2£ 2 ) ytanh(y) + 7sinh 2 (?/) 



One can see that again the perturbation matrices contain elements of different order in e and 
therefore the calculation of energy levels needs first to derive the perturbed eigenfunctions 
to a relevant order in e as it has been done in the Appendix C. 

To the relevant order in e, the eigenvalue equations for the scalar wave function com- 
ponents sj^(y), s7(y) follow from the matrix elements of the operator in Eq. (|D.l| ), 



(4 - E s ) + [4e 2 tanh 2 (y) + A} 1 

+4e (1 + f) sech 2 (y) sinh(y) sj(y) = 
(d+ - E s ) sj(y) + [4e 2 sech 2 (y) + A 22 



sj(y) 
sj(y) 



+4e (1 + £) sech (y) sinh(y) s1j~(y) 







(D.4) 



These equations have terms of different order in e in two sectors of solutions. 

Actually, for even-odd solutions with s^(y) = (j)f(y), the upper even component 



-(0) 



+ e 



+ C(e 4 ) contains the zero mode ( |5.10| ) of the operator d$ as a leading 
term. Meanwhile, the lower odd component <p~ has an order of e and coincides with the 



unperturbed one - see Eq. ( J5.10 ) in its functional form. Indeed, as one expects that the 



lightest eigenvalues are of the order E s 
Eq. flDl]) reads, 



~ e 



then the dominating part of the second 



d 



--(0) 



+ 4e (1 + f ) tanh(y) sech(y) </> +(0) = 0. 



(D.5) 
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Comparing this equation with the unperturbed one, corresponding to £ = 0, from Eq. ( 5. lCj ) 
one immediately finds that 



~ (0) (2/) = ~\ ~T e (1 +0 sinh(y)sech 2 (y) 



(y) 



(D.6) 



On the other hand, as d 

-(i) 



-(0) 



the first Eq. (|D.4| ) determines the next-to-leading 



contribution </> ■ into the even part of the wave function. Now this equation allows us to 
obtain both the approximate eigenvalue and the correction term </>j • After projection 



onto the zero-mode 



-(0) 



of the operator cQ~, one gets rid of the unknown function <j)- 
and calculates eventually the branon mass: namely, 



(i) 



2 r+°° 
m{) = /3 2 ^ = /3 2 / dy 

J — oo 



4e 2 tanh 2 (y) + A} 1 



-(0) 



(v) 



+4e (1+C) sech 2 (y) sinh(y) 



(2/)0 +(O) (y)} 



(D.7) 



In the sector of odd-even solutions (y) = hj (y) the lightest-state wave function has the 
opposite order in e, i.e., 

hj i0 \y) ~ e, h+(y) = + e 2 h +(i) + 0(e 4 } 

Thus the role of two Eq.s ( p.4| ) is merely interchanged. In particular, it can be easily 
derived that 



- e(l + Oysech 2 (y)/V2 = h-^(y) 
h + (°\y) = sech(y)/v / 2, 



(D. 



to be compared with Eq. ( |5.12[ ). Respectively, when we project onto the zero-mode hj (y) 
of the operator d^, the second equation helps us to obtain the Higgs boson mass, 



in 



2 r+oo 

3 \ _ al t? «2 



dy 



4e 2 sech 2 (y) +A 22 1 \h +( -°\y) 



+4e(l+0 sech 2 (y)sinh(y) /i~ (0) (y) 



fD.{ 
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